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Abstract 

We have developed a Mathematica package capable of performing T-matrix algebra 
in arbitrary (integer) dimensions. As an application we can compute Fierz transfor- 
mations. 
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COPYRIGHT AND DISCLAIMER OF WARRANTY 



Copyright © 2001 Ulf Gran. 

This documentation is part of the GAMMA package, available for download at 
http://fy.chalmers.se/~gran/, developed for Mathematica version 4.0. 

The GAMMA package is distributed in the hope that it will be useful but WITHOUT 
ANY WARRANTY. Neither the author nor any distributor accepts responsibility 
to anyone for the consequences of using it or for whether it serves any particular 
purpose or works at all, unless he says so in writing. 

Permission is granted for use and non-profit distribution of this package, in its origi- 
nal form, providing that this notice be clearly maintained, but the right to distribute 
this package for profit or as part of any commercial product is specifically reserved 
for the author. The same permission and reservation applies to the formatted doc- 
ument distributed with this package. 

If you find the GAMMA package useful for doing research, acknowledging the au- 
thor's work by proper reference to the package would be appreciated. 
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1 Introduction 



The motivation to develop this package arose from calculations in eleven-dimensional 
supergravity p], 0], which required an immense amount of both pure T-matrix al- 
gebra and Fierz transformations to be performed. The rule-based programming of 
Mathematica || is ideally suited for handling algebraic computations, compared to 
ordinary procedural programming. For a very nice comparison between rule-based 
and procedural programming, as well as a toy implementation of T-matrix algebra 
in four dimensions, see [Q], which describes a Mathematica package performing T- 
matrix algebra in (4— e) dimensions. This was the first Mathematica implementation 
of T-matrix algebra according to the 't Hooft-Veltman scheme. There are quite a few 
programs handling T-matrices in four-dimensions []5|, ^, [7|, ||, ||, 1C, 11], intended to 



be used when computing Feynman diagrams. There is also a Mathematica package 
handling computations in D = 10, N = 1 supergravity, including T-matrix algebra, 

In section ^, we will introduce and explain the functions in the package concerning 
the algebraic manipulation of T-matrices. This will be done through simple exam- 
ples. In section || we will turn to the functions concerning fierzing. A complete list of 
the functions in the package, including brief descriptions, can be found in appendix 
|A|. In appendix P we list the built-in Mathematica functions that are redefined by the 
package. GAMMA is available for download at http://fy.chalmers.se/~gran/. 



2 T-matrix algebra 

To load the package type 



In[l] :=«GAMMA.m 



and remember to include the path if the package is put in a directory where Math- 
ematica does not automatically search. By default, the space-time dimension is set 
to 11 and the spinor dimension is set to 32. If you want to use any other values they 
must be specified before any calculations can be done. This is done through, e.g., 



In[l] :=SetDim[10] ; SetSpinorDim [16] ; 



The spinor dimension is used only when calculating the trace of T-matrices. We now 
introduce the various objects that can be handled. Kronecker deltas, T-matrices, 
tensors and tensor-spinors are entered like 



In[l] :=Delta[{a,b},{c,d}] 
0ut[l] 
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In [2] :=GammaProd[{a,b},{c}] 

out [2] =r ab **r c 

In [3] : =Tensor [X , {a , b , c}] 
Out [3] =X abc 

In [4] : =TensorSpinor [Y , {a , b , c}] 
Out [4] =Y abc 



The indices in each list in these expressions are assumed to be antisymmetric and will 
be put in a canonical order (defined by Mathematica's sort algorithm, which there- 
fore means alphabetical order), sometimes resulting in a sign in front of the canon- 
ically ordered expression. In order to discern tensors from tensor-spinors (since we 
do not use explicit spinor indices) we have added a tilde over the symbol designating 
the tensor-spinor. This indicator can be changed by, e.g., 



In [1] : =SetTensorSpinorInd [OverHat] 



in which case the Mathematica function OverHat is applied to the symbol designating 
the tensor-spinor. See the Mathematica documentation for a complete list of possible 
choices. It is also possible to change the formatting of the symbols designating a 
tensor, where the default function is 



In[l] :=SetTensorInd [Identity] 



Since all calculations are assumed to be performed in a Lorentz frame, we do not 
have to care about whether an index is upstairs or downstairs and contracted indices 
may both be upstairs or downstairs. Enforcing explicit antisymmetry in some indices 
in an expression like T[ a X b ] is done by 



In[l] :=ASym[GammaProd[{a}] Tensor[X,{b}] ,{a,b}] 

Out[i]=|(x b r a -x a r b ) 

I I 

Using the function Sym[] instead enforces symmetry in the specified indices. 

We now turn to the functions concerning the manipulation of T-matrices. We 
have already seen how to write a product of T-matrices, if we want to expand a 
general product we write 



In [1] : =GammaExpand [GammaProd [{a,b} , {c} , {d}] ] 

out [i]=#?r a **r d - £*r b **r d + r abc **r d 
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where the product of the first two T-matrices have been expanded. Since only the 
product of two T-matrices are expanded each time GammaExpand [] is applied, re- 
peated application is sometimes needed. In the expansion we only use the algebra 
{r a ,r;,} = 25% and therefore, as long as the Kronecker deltas have an equal num- 
ber of upstairs and downstairs indices, there will be no difference between various 
signatures. In a specific expression, if we want a Kronecker delta to have, e.g., only 
upstairs (or downstairs) indices we get a trivial difference since using Minkowski 
signature the Kronecker delta becomes an rj while it stays a delta in the euclidean 
case. Therefore, with the understanding that a DeltaProd [] has an equal number of 
upstairs and downstairs indices, the expressions will look the same for all signatures. 

Since we know that the expression we started with in this example above is 
antisymmetric in a and b we can simplify the result by putting these indices in 
canonical order by 



In [2] : =ACanonicalOrder [% , {a , b}] 

out [2] =-2^r b **r d + r abc **r d 



remembering that that the expression is antisymmetric in a and b. Even though the 
simplification in this simple case is not that significant the importance to use this 
function, when applicable, in more complicated cases can not be underestimated. 
The function SCanonicalOrder [] puts indices in canonical order assuming symme- 
try instead. When indices are contracted between adjacent T-matrices it is more 
efficient to use 



In[l] : =GammaContract [GammaProd[{a,b,c},{b,c}]] 
Dut[l]=-90r a 



which is specialized for this situation, than to use GammaExpand [] . To compute 
traces of expressions involving T-matrices we write0 



In [1] : =GammaTrace [GammaProd [{a, c} , {b, c}] ] 
Out [1] =-320<5§ 



When using tensor-spinors, which are assumed to be irreducible and therefore T- 
traceless, it is sometimes needed to rearrange the T-matrices in order to be able to 
use the T-tracelessness to simplify an expression. We see in this example 

2 When tracing a T-matrix with the same number of indices as the spacetime dimension, when 
this number is odd, GammaProd [] is kept and the user can himself replace it with a Levi-Civita 
tensor density according to his conventions (note that GammaProd [] is not multiplied by the spinor 
dimension) . 
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In [1] : =GammaProd [{a , b , c}] **TensorSpinor [Y , {c}] 
Out[l]=r abc **Y c 



In [2] : =GammaExtract , {a , b}] 
Out [2] =-2r a **Y b 



that we have to use GammaExtract [] to be able to simplify the expression. What 
GammaExtract [] does is to extract one T-matrix at a time, having an index con- 
tracted to the tensor-spinor, from, in this case, the r abc term, in this way enabling 
the use of T-tracelessness, e.g., T C Y C = 0. Note that it is very important to use 
**, i.e., NonCommutativeMultiply [] , in this example since there are more than one 
object carrying spinor indices. After some manipulations, one often ends up with 
expressions where contracted, i.e., dummy, indices are denoted by different symbols 
in otherwize identical expressions like in 

i 1 

out[i]=x a r a + x b r b 

i i 

In this case Mathematica can not tell that these terms are identical and we have to 
rename the dummy indices using some canonical set of indices. The is done by 



In [2] : =RenameDummy [%] 
Out [2] =2X dl r dl 



The set of indices {dl, d2, . . .} are used to denote dummy indices and are therefore 
preferably avoided. It is possible to cut, paste and copy formatted expressions 
containing T-matrices and also to enter expressions as they look in the output form. 
Converting expressions to TgX format is done by TeXForm [expr] . 

3 Fierz transformations 

We can also use GAMMA to perform Fierz transformations. This will be illus- 
trated using a simple example in which the free vector indices of the T-matrices are 
contracted to a tensor. Start by making the definitions 

i 1 

structures={xl , x2 , yl , y2 , y3 , zl , z2 , z3} ; 

ind=HoldForm[|tr[P**X]**Q + §P**X**Q] ; 

xl={{a},{b},ind,{{},{}}}; 
x2={{cl},{cl,a,b},ind, {{},{}}}; 
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yl={{a,b},{},ind, {{},{}}}; 
y2={{a,cl},{b,cl}, ind, {{},{}}}; 
y3={{cl,c2}, {a, b,cl,c2},ind, {{},{}}}; 

Zl={{a,b,cl,c2,c3},{cl,c2,c3},ind, {{},{}}}; 
z2={{a, cl , c2 , c3 , c4} , {b , cl , c2 , c3 , c4} , ind, {{} ,{}}}; 
z3={{cl , c2 , c3 , c4 , c5} , {a, b , cl , c2 , c3 , c4 , c5} , ind ,{{},{}}}; 

i i 

The list structures contains all the different terms we want to relate by fierzing. 
A term wl={{a} , {b} , ind, {{c} , {d}}} corresponds to 

(T%,(TV§. (1) 

Note that the terms in this example have two free vector indices, a and b. To perform 
the Fierz transformation we contract two spinor indices of each term with a basis 
of T-matrices, and for symmetric spinor indices and C, T^ 3 ' and T^' 

for antisymmetric spinor indices. In our case we are interested in three symmetric 
indices, (afirfS. Since we do use explicit spinor indices on the T-matrices we use the 
expression ind to specify the symmetry of the four spinor indices. In ind, P will be 
replaced by the first T-matrix in a term and Q by the second. X will be replaced by 
the basis matrices and tr by the function GammaTrace. It is then easy to see that 
ind corresponds to the symmetry specified above^]. We are now ready to perform 
the Fierz transformation, 



In [1] : =FierzSolve [Fierz [{yl ,y2 ,y3} , structures , {a,b} , S , 
MakeTensor] ] 

Out[l] = 
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where the first argument to Fierz is the terms that are chosen as the basis in 
which the other terms will be expressed. It is important to know exactly how many 
linearly independent terms there are and to get this information we recommend the 
lie algebra program Lie O]. The second argument is the list of terms, the third is a 
list of the free vector indices, the fourth is S since we are fierzing symmetric indices 

3 Contract a basis T-matrix, T x , with two of the symmetrized indices, (T x ) al3 T^ al3 r^ g = 
|tr-(r x r p )r^5 + |(r p r x r Q ) 75 , where the RHS is the expression to be used as ind. 
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in this example (use A for antisymmetric indices). The last argument MakeTensor 
is used to indicate that the symmetry of the free indices is that of a tensor. In the 
output, each row contains the coefficients for a term in the chosen basis. In this 



case, e.g., z2 = -f yl + §y2 + §y3. 



When we Fierz tensor-spinors, care has to be taken in order to be able to use 
the T-tracelessness. If, e.g., the index b in the previous example is contracted to a 
vector-spinor having the spinor index 7 we would have to use 



ind=HoldForm [|tr[P**X]**Q**TensorSpinor[Y, {b}] 
+ §P**X**Q**TensorSpinor[Y, {b}]] 

1 1 

The tensor-spinor is automatically removed after the simplification is complete when 
fierzing. 

Fierz transformations like the one above where the free vector indices are con- 
tracted to an antisymmetric tensor are relatively simple to perform, it is harder to 
Fierz hook-tensors^. Even though GAMMA may not be able to directly calculate 
the Fierz transformation as shown above, it can still be of great use. We can try to 
use Fierz [] but if that fails we can use FierzTerm[] instead, which contracts the 
basis of T-matrices to a specific term, or list of terms, simplifies and then returns 
the result. The problem with hook-tensors is that a lot of Kronecker delta func- 
tions are generated and it is not easy to rewrite them on a canonical form allowing 
Mathematica to determine which terms are linearly independent. 

The way to specify a hook-tensor term is a bit different from the tensor case. 
Let us look at a term arising when we want to Fierz a (11000) representation, 

(T c ^) a(3 (Tf a )ysZ ClC2 , a , (2) 

where we have one free vector index. This term is represented by 



f 1={{} , {f },{cl,c2}, {a}, 2, ind, {{},{}}}; 

1 1 

Here, the first entry corresponds to the vector indices on the first T-matrix, not 
including indices contracted to the hook-tensor, and the second entry is defined in the 
same way but for the second T-matrix. The third entry is the antisymmetric hook- 
tensor indices, except those that sit on a Kronecker delta function (after removing the 
hook-tensor^]), and the fourth entry is the symmetric index on the hook-tensor. The 
fifth entry is the number of vector indices on the first T-matrix, ind again specifies 
the structure of the spinor indices and the last entry is the potential Kronecker 
delta function indices. When we remove the hook-tensor in (|J) we must enforce 
the hook-tensor symmetries on the T-matrix terms, this is done by the function 

4 A hook-tensor is denoted Z cl ... Cp ,a satisfying Zr ci .„ c tt i = and Z cl ... Cn _ lc c = 0. 
5 For a term as in (||) we have to introduce a Kronecker delta function when removing the 
hook-tensor if the free index f would sit on the hook-tensor instead of on a F-matrix. 
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MakeHook [term.gmatrix] , which at the same time contracts two spinor indices to 
the spinor indices on the T-matrix gmatrix as specified by ind. By writing 



In [1] : =FierzTerm [{f 1} , {f } , S , MakeHook] 

i i 

we generate the expressions obtained by contracting the symmetric T-matrix basis, 
in this case, to the hook structure specified by f 1, treating the indices contained in 
the list of the second argument as antisymmetric. These expressions can then be 
studied "by hand" and in this way we can also simplify hook fierzes. 



4 Summary 

We have presented a Mathematica package capable of performing T-matrix alge- 
bra in arbitrary (integer) dimensions. The various functions have been introduced 
through examples and a complete list of all functions in the package, including brief 
descriptions, can be found in appendix As a specific application, we have consid- 
ered Fierz transformations, which in general require a lot of T- matrix calculations. 
Note that some of the built-in functions in Mathematica have been altered, see 
appendix [B| for a complete list. 

Comments, suggestions for improvements and reports regarding any discovered 
bugs are appreciated. When reporting bugs, please include a Mathematica note- 
book illustrating the bug in question. Reporting bugs in this way will result in the 
quickest possible bug-fixes. Please note that only the most basic error handling is 
implemented. Improved error handling invariably results in a slower program and in 
the choice between error handling or speed, we chose the latter. This unfortunately 
means that what might at first sight look like a bug quite often can be attributed 
to some error in the input. 
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A List of commands 

Here is a list of all the user-available functions in the GAMMA package together 
with a briefs descriptions. 

ACanonicalOrder [expr, list] assumes that expr is antisymmetric in the indices 
given by list and puts them in canonical order, thus enabling additional 
simplification. 

ASym [expr, list] explicitly antisymmetrizes expr in the indices contained in list. 
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Delta [list l,list2] represents a Kronecker delta with the two groups of indices 
given by listl and list2. 

Fierz [basislist , structures , asymlist , sym, function] performs the Fierz trans- 
formation of the terms specified in structures using the terms in the list 
basislist as basis elements. Antisymmetry is assumed in the indices con- 
tained in the list asymlist and sym (A or S) gives the symmetry of the spinor 
indices to be fierzed. Finally, function is applied to generate explicit T-matrix 
expressions given a term in structures. 

FierzSolve [eqlist] solves the list of equations generated by Fierz [] . 

FierzTerm [termlist , asymlist, sym, function] as Fierz [] but does not perform 
the complete Fierz transformation. Instead it produces a list of expressions 
obtained by contracting two of the spinor indices, as specified by ind, with the 
basis of T-matrices specified by sym. 

GammaContract [expr] simplifies the expression expr containing contracted indices 
on adjacent T-matrices. 

GammaExpand [expr] expands the product of the first two T-matrices in expr. 

GammaExtract [expr, list] extracts T-matrices that contract indices on a tensor- 
spinor, assuming antisymmetry in the indices contained in list, in order to 
be able to invoke the T-tracelessness of the tensor-spinor. 

GammaProd [listl , . . . ,listn] represents the product of n T-matrices with indices 
given by listl to listn. 

GammaTrace [expr, list] computes the T-trace of an expression expr containing 
T-matrices, enforcing antisymmetry in the indices contained in list. 

MakeHook [term,basiselement] makes an explicit T-matrix expression from term 
contracting the T-matrix basiselement with two of the spinor indices, as 
specified by ind. 

MakeTensor [term, basiselement] makes an explicit T-matrix expression from term 
contracting the T-matrix basiselement with two of the spinor indices, as spec- 
ified by ind. 

RenameDummy [expr] renames dummy (contracted) indices in expr, thus enabling 
additional simplification. 

SCanonicalOrder [expr , list] assumes that expr is symmetric in the indices given 
by list and puts them in canonical order, thus enabling additional simplifi- 
cation. 

SetDim [dim] sets the space-time dimension to dim. 

SetTensorlnd [function] sets the Mathematica function to be applied to the sym- 
bol designating a tensor to function. 
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SetTensorSpinorlnd [function] sets the Mathematica function to be applied to 
the symbol designating a tensor-spinor, in order to discern it from a tensor, to 
function. 

SetSpinorDim [dim] sets the spinor dimension to dim. 

Sym[expr,list] explicitly symmetrizes expr in the indices contained in list. 

Tensor [sym, list] represents a tensor, denoted by the symbol sym, with indices 
given by list. 

TensorSpinor [sym, list] represents a tensor-spinor, denoted by the symbol sym, 
with (vector) indices given by list. 

B Redefined Mathematica built-in functions 

NonCommutativeMultiply [] : Redefining NonCommutativeMultiply [] is what en- 
ables us to use new, noncommutative, objects like T-matrixes and spinors in 
Mathematica. NonCommutativeMultiply [] can also be entered as **. 

Power [] : In order to be able to interpret expressions entered in the output format, 
and also to be able to cut, copy and paste expressions, we have to modify 
Power [] so that e.g., P 1 is replaced by GammaProd [//] . 

Subscript [] : The same reason as for Power [] except that Subscript [] takes care 
of downstairs indices instead of upstair ones. 
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